(N 



in 



> 



Inter-particle potentials in a scalar QFT with a Higgs-like 

mediating field 

Alexander Chigodaev and Jurij W. Darewych 
Department of Physics and Astronomy 
York University, Toronto, ON, M3J 1P3 



O ■ January 20, 2013 

, Abstract 



We study the inter-particle potentials for few-particle systems in a scalar theory with a non-linear 
mediating field of the Higgs type. We use the variational method, in a reformulated Hamiltonian formal- 
ism of QFT, to derive relativistic three and four particle wave equations for stationary states of these 
systems. We show that the cubic and quartic non-linear terms modify the attractive Yukawa potentials 
but do not change the attractive nature of the interaction if the mediating fields are massive. 



P-i- Introduction 



Inter-particle interactions in systems described by QFTs, such as QED or the Yukawa model and its 
variants (like the boson exchange model in nuclear theory), are specified by the interaction term (or terms) 
\l . of the underlying Lagrangian density. For theories where the mediating fields are linear (e.g. the Electromag- 

netic field for QED or the boson fields in nuclear theory), the relevant interaction terms of the Lagrangian 
(density) are of the generic form ^^A^, where is the particle field (charged fermions in QED; nucleons 
' in nuclear theory) and A describes the mediating field (photons in QED and mesons in nuclear theory). 

. Few particle wave equations that stem from theories with such interaction terms in the Lagrangian density 

.• ■ ^ _ ... . . . e"^*" 

' are, in the non relativistic limit, Schrddinger equations with inter-particle potentials given by ±a , 

where a is the coupling constant, /i is the mass of the mediating-field quantum and r is the inter-particle 
distance. However, there are theories in which the mediating fields are non-linear, such as, for example, in 
the Higgs model, QCD and other models. In these theories, the interaction part of the Lagrangian density 
includes terms that are typically of degree 3 and 4 in the mediating field A. The analytic derivation of the 
. inter-particle potentials in such cases becomes a more difficult and generally not achievable task. 

In this paper, we shall address this problem for a scalar field theory with a non-linear mediating field, 
specified by the Lagrangian density (units: h = c = 1): 

C - d'-r d,^ -m'ct)* ct^-gxct>*ct>-\ {^*4>f + \ {d'^x d,x ~ X^) ~ \ri ~ \o ■ (1) 

The parameters m and ^ represent the bare masses of the scalar and mediating field quantum respectively, 
rj are coupling constants with dimensions of mass, while cr, A > are dimensionless coupling constants. 
Note that the x field is of the form of the Higgs field of the Standard Model. Moreover, the non-linear terms 
in the field x are of a form that mimics those of the gluon field in the QCD Lagrangian ii fi = 0. The A term 
in the Lagrangian density is required to keep the classical ground energy bounded from below. It leads to 
a repulsive contact (delta function) inter-particle interaction for few-particle systems [1,2]. Henceforth, we 
set A = since it has negligible effect on what follows in this work 

The aim in this paper is to derive relativistic wave equations for few-particle states of this model, with an 
emphasis on the study of the inter-particle interaction potentials. We shall use a reformulated Hamiltonian 
formalism, together with the variational method, to accomplish this task. This approach has been used 
effectively in the study of few-boson states in the scalar Yukawa theory [3, 4] and relativistic few-fermion 
bound states in relativistic quantum mechanics [5] and QED [6,7]. 
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Reformulation 



Wc introduce a reformulation of the Lagrangian density with the intention of expressing the effect of 
the mediating field x through its propagator. The reformulation corresponds to a partial solution of the 
equations of motion. The equations of motion corresponding to the Lagrangian (1), derived from the action 

principle 5 J dx C = 0, are 

{d^ + n?) (j>{x) = -g(p{x) x{x), (2) 
-g^*{x)cP{x)-rjxHx)-<TxHx). (3) 

Equation (3) has the integral representation {i.e. "formal solution") 



X{x) = xo(x) + J dx' D{x - x') p{x'), (4) 

where x = {t,x), dx = dtdx, and 

p{x) = -g<i)* {x)(t){x) -r]x^ (x) - ax^ {x) (5) 

is the "source term" of the inhomogeneous equation. The function xo{x) satisfies the homogeneous (or free 
field) equation with the right hand side of (3) equal to zero, while D{x — x') is a covariant Green function 
(or the Feynman propagator) for the mediating field X; such that 

{d'^ + fj."^) D{x - x') = d^^\x - x'). (6) 

Recall that the Green function can be expressed as 

Dix - X') = Dix,x') = lim/ ^.^^^^r''''-"'- 

We shall not be concerned with processes that involve free quanta of the mediating field XO) and so 
henceforth we shall leave it out. Then, substituting (4) into (2), omitting terms containing xoi we obtain 
the equation 

[d'^+m^) (pix) = g^(l>{x) j dx'(l>*{x')(l){x')D{x,x') 

+ gil(t^{x) j dx'x^{x')D{x,x')+ga^{x) j dx' x\x') D{x,x'). (8) 
Equation (8) is derivable from the Lagrangian density 

C = d^(l>*{x) d>'(t>{x) - (t>* {x)(l){x) + \9^ j dx'4>*{x)(t>{x)D{x,x')(t)*{x')<t>{x') 

+ gr,r{x)cp{x) J dx'D{x,x')x\x')+gar{x)cj>{x) J dx' D{x,x')x\x') (9) 

provided that the Green function is symmetric, i.e. D{x — x') = D{x' — x). 

For the case of a linear mediating field, i.e. when r] = a = 0, the reformulated Lagrangian density (9) 
gives field equations that involve the particle fields only; the mediating field x appears only through the 
propagator D{x,x'). The reformulated Lagrangian density (9), with = cr = 0, is convenient for the study 
of few-boson relativistic bound states in the scalar Yukawa theory as will be pointed out below and as was 
shown in references [3,4]. 

However, for the present case of a non-linear mediating field, i.e. r],a ^ 0, it is not possible to obtain a 
closed-form solution of (4) for the field x in terms of D{x, x') and (j), thus we need to resort to approximation 

schemes. An obvious one is an iterative sequence, based on (4). The first order approximation corresponds 
to substituting the formal solution (4), with p = —g(j)*(p (with xo left out), into the Lagrangian density (9). 
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This yields the first-order iterative approximate expression 

L = d^(j)*{x) d^(l){x) - (p* {x)(j){x) + ^9^ J dx'(j)*{x)(l){x)D{x,x')(l)*{x')(j){x') 

+ gr](l)*{x)4>{x) j dx' dz' dz" D{x,x') D{x' , z') D{x' , z") p{z') p{z") 

+ ga(l}*{x)(i){x) j dx'dz'dz"dz"' D{x,x')D{x',z')D{x',z")D{x',z"')p{z')p{z")p{z"'). (10) 

The interaction terms in the Lagrangian density (10) do not contain the mediating field x explicitly but 
rather implicitly through the mediating field propagators. The advantage of reformulating the Lagrangian 
density to the present form will be evident when we implement quantization. As we shall sec, this will 
enable us to use simple few-particle trial states free of the x field quanta while still probing the effects of all 
interaction terms, including the non-linear terms {i.e. those with r],a ^0). 

The Hamiltonian and momentum densities corresponding to (10) are obtained from the energy-momentum 
tensor in the usual way: 

T'^" = ^7Tr—.d''<i>i-9^" (11) 
a(^^^<pi) 

where the index i = 1,2 stands for the fields (f>i = (j> and (j)2 = 4>* ■ The component of the energy- 
momentum tensor is the reformulated Hamiltonian density: 

?^ = </)n^ + n^. - £ = + + Ui, + Hi„ (12) 

where 

= n^. + ( v<^* ) • ( v<^) + (13) 
1 



^/i = - ^ 9^^*{^)<t>{-r) j dx'r{x'mx')D{x - x'), (14) 
= -gn<f)*{x)<j){x) J dx'dz'dz" D{x,x')D{x',z')D{x',z") p{z')p{z"), (15) 
Hi, = -ga(j)*{x)(j){x) J dx' dz' dz" dz'" D{x,x') D{x' , z') D{x' , z") D{x' , z'") p{z') p{z") p{z"'). (16) 

dC ■ dC ■ 

The conjugate momenta of the fields are defined in the usual way, i.e. 11^ = — r = 4>* and 11^. = — = (f). 

d(j) dip* 

Observe that every term of the reformulated Hamiltonian density has dimensions of (as it should). Notice 
that the Hamiltonian density is non-local, i.e. it depends on more than one space-time variable. Equation 
(12) is the Hamiltonian density that we shall use in the present work. 

The T*^' components of the energy-momentum tensor are the momentum density components 

V' = -n^diif>-n^.di(p\ (17) 

Note that the interaction terms do not contribute to the momentum density above. We shall make use of 
this expression when the total momentum of the few-particle systems will be discussed. 



Formalism and Quantization 

We will work in the Hamiltonian formalism of QFT where the basic equation to be solved is the 4- 
momentum eigenvalue equation 

P'' 1^-) = l^r), (18) 

where P" is the energy- momentum operator of the quantized theory which follows from equations (12) (as 
specified below) and (17), Q"^ = {E,Q) is the energy-momentum eigenvalue and \^) are the corresponding 
eigcnfimctions. It is not possible to obtain exact solutions for the v = component of equation (18). 
Consequently, if one wished to solve for the energy of a system in question, then approximation methods 
must be used. For u = 1, 2, 3 the solutions of (18) are free-field-like as we point out below. 
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Variational approximations to the i/ = component of the eigenvalue equation (18), i.e. the Hamiltonian 
component, require the evaluation of 

S{%\H-E\^t)=0, (19) 

where arc suitably chosen trial states that contain adjustable features (functions, parameters). The 
applicability, accuracy and usefulness of a variational approximation depend on the choice of the trial states. 
They must be appropriately chosen for every case being studied. 

We now proceed to the canonical quantization of the theory with the Hamiltonian density (12). First, 
we expand the fields in terms of its Fourier modes: 

c/,{x) =Jdp [(277)3 2a;p] (a (p) 6"^^- + 6+ (p) e'f-) , (20) 

r{x)= ydp[(27r)3 2a;p]"^(a^(p)e^f-" + 6(p)e-*f"), (21) 

where cjp = + p^ . Note that since the mediating field x docs not appear in the Hamiltonian explicitly, 
we need not be concerned with its Fourier mode representation. The next step is to promote the fields to 
the status of operators and impose the following non- vanishing, equal time commutation relations: 

[^{x),U^{y)] = [0*(.xO,n<j.(y)] =2(5(x-y). (22) 

All other commutators of the field and conjugate momentum operators vanish. The operators a(p) and 
6(p) are the free particle and anti-particle annihilation operators respectively, whereas the operators a\p) 
and b^{p) are the free particle and anti-particle creation operators. The vacuum state |0) is defined by 
a(p)|0) = &(p)|0) = 0. Note that the canonical quantization is performed in the interaction picture. 

Expressing the creation and annihilation operators in terms of the fields and its derivatives, the commu- 
tation rules become 

[a(p), at(q)] = [6(p), b\ci)] - S{p - q). (23) 

In the present work, we shall not be concerned with vacuum energy questions so we shall normal-order 
the creation and annihilation operators in the Hamiltonian. To obtain the Hamiltonian operator we integrate 
out the spatial coordinates from the Hamiltonian density (12): 

H{t) = j dyi: n{t,x) :, (24) 

and express it in terms of the particle and anti-particle operators a, b, and . This is straightforward 
for the free-field part of the Hamiltonian as can be seen in standard books on QFT. However, obtaining the 
analogous expressions for the interaction terms is tedious and this has been accomplished with the use of a 
computer code. 

The momentum operator is similarly obtained by normal-ordering and integrating out the spatial coor- 
dinates from the momentum density (17): 

dx : V{x) : . (25) 

Notice that since the momentum density operator does not involve the interaction terms, it is time indepen- 
dent. 

For the description of stationary few-particle bound states, it is convenient to switch to the Schrodinger 
picture where the time-dependence of the Hamiltonian is removed. The two pictures are related by 

|*/(t))=e^^«*|*s(i)) (26) 
where Hq is the free-field part of the Hamiltonian, i. e. with interactions turned off. 

Particle-Antiparticle State 

The simplest particle-antiparticle trial state can be expressed in terms of Fock-states as follows: 

|^E'2) = J cipi,2 i^(Pi,2) a+(Pi) b\p2) |0), (27) 
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where -F(pi,2)) is an adjustable coefficient function which is determined variationally. The subscript notation 
means that dpi,2 = c?Pi dp2 and F{pi,2) = F(pi,p2); it is used throTighout this paper. 

To implement the variational principle (19), the matrix element {'^2\ H — E \^2) is worked out and is 
varied with respect to the adjustable coefficient function F* . This leads to the following relativistic wave 
equation in momentum space describing a particle-antiparticle system: 

F{P1,2) Ki +'^V2-E)= j dp[^2 3^2,2 (Pi,2, Pl,2) i^(p'l,2), (28) 

where 3^2,2 is the relativistic Yukawa particle-antiparticle interaction kernel (inter-particle "momentum- 
space" potential). In the interaction picture, it is given by 



^2,2(pi,2, Pl,2) = o /Lx3 /", " ,,—r- '^(P'l + P2 - Pi - P2) 

8(27r)-^ yWp/Wp/Wp.Wp^ 

f 1 1 



2 '("p'i+"p^-'^Pl-'^P2) 



ip[ -Pl)2 ^ - ip'2-P2r ^ - {Pl +P2y ^ - {p[ +p'2ff' ^^^^ 

Notice that the first two terms in the square brackets on the right hand side of equation (29) correspond 
to onc!-quantum exchange and the last two to virtual annihilation Feynman diagrams. These diagrams are 
shown in Figure 1. 

The time dependence of the particle-antiparticle interaction kernel 3^2,2 can be "rotated away" by the use 
of equation (26). If we start with the trial state j^l/J,) = e'^''*|vl>2) then the corresponding interaction kernel 
in the Schrodinger picture will be time independent, i.e. we obtain equation (29) with t = 0. Therefore, the 
Schrodinger picture (time independent) interaction kernels can be found, in effect, by setting t = in all 
interaction picture kernels. This is what we shall do in the remainder of the paper. 

Unfortunately, the quantized version of the non-linear terms (15) and (16) of the Hamiltonian density are 
not probed by the simple particle-antiparticle trial state (27), i.e. the matrix elements (^2! Hj^ \^2) vanish 
for i = 2,3. To examine the effects of these terms we must use more elaborate trial states than (27). As it 
is, with the trial state (27), the problem reduces to the scalar Yukawa model which has been discussed by 
Ding and Darewych [3]. Nevertheless, it is useful to recount some details. 

For the particle-antiparticle trial state (27) to be an eigenstate of the momentum operator (25) we require 
that P 1*2) = Q 1*2)- In the rest frame, Q = 0, this requires that i^(pi_2) = (5(pi + P2) /(Pi)- The centre 
of mass motion separates and the relativistic momentum-space particle-antiparticle wave equation simplifies 
to 

/(P)(2^p-^)-^|^/^/(p') 

' ' ' ^ +U-^±-A\. (30) 



[ + (p' - p)2 - (wp, - a;p)2 2Vm'-4w2; 2\^/x2-4a;2^ 

The first term on the right hand side of equation (30) corresponds to the one-"chion" exchange {i.e. one me- 
diating field quantum exchange) while the second corresponds to the virtual annihilation Feynman diagrams. 

In the non-relativistic limit (i.e. p^ << m^) the Fourier transform of equation (30) yields the expected 
Schrodinger equation in coordinate space for the relative motion of the particle-antiparticle system: 

--V2V(x) + y(a;)V(x) =eV(x), (31) 

m 

where e = E — 2m and the potential energy depends on the particle-antiparticle distance: 

p— 4 TT 

n-) = -.^ + ;^^(x), (32) 

9^ 

where a„ = — 7; is a dimensionless coupling constant. The first term, due to the one-chion exchange, 

16 7rm^ 

is the usual Yukawa potential and is always attractive {i.e. gravity- like) in this scalar theory. The second 



5 



T 




(a) (b) 



Figure 1: The one-chion exchange (a) and virtual annihilation (b) diagrams of the Yukawa kernel equation 
(29). 

term, due to the virtual annihilation, can be either attractive or repulsive depending on the values of the 
parameters fj, and m. It is a correction to the Yukawa potential and is a feature of the quantum field theory. 

The relativistic equation (30) is not analytically solvable, so approximation methods must be used. 
Perturbative and variational approximations are presented in the paper by Ding and Darewych [3], along 
with a comparison to results obtained using the ladder Bethe-Salpeter equation and various quasi-potential 
approximations. 



Three Particle State 

To observe the effects of the non- linear terms of the Hamiltonian density (12) on the inter-particle 
potential, wc must consider trial states with either more particle content or more Fock-space components. 
We carry out the easier task first and consider a three identical particle trial state given by 

1*3) = j dpi,2,3 i^(pi,2,3) a^(pi) at(p2) at(p3) |0), (33) 

where F(pi 2,3) is a three-particle function to be determined variationally. Note that this trial state 
can be taken to be an eigenstate of the momentum operator (25), i.e. P 1^3) = Q |^'2) with the choice 
-F'(Pi.2..3) = S{pi + P2 + P3 — Q) /(pi,2): where Q is the constant total momentum of the three particle 
system. Consequently, the wavcfunction will be of the form where the centre of mass motion is completely 
separable, just as we saw for the particle-antiparticle case. 

The addition of an extra third particle operator to the trial state significantly complicates the evaluation 
of the matrix element. However, the implicit symmetry of equation (33) with respect to interchanges of 
momentum variables due to the identity of the particles makes the task easier. It permits us to carry out 
the evaluation of the matrix element for the three particle trial state in terms of the symmetrized function 
Fs: 

6 

i^s(pi,2,3)= Yl PiPn,i2,is), (34) 

where the summation is on the six permutations of the indices 1, 2 and 3. Making use of this symmetrization 
enables us to keep the arguments of the kernels non-permuted while storing all information about the 
symmetry under interchanges of the momentimi variables in Fg. 

We then calculate the matrix element (^3! H — E l^'s), carry out the variational derivative with respect 
to F* and set it to zero (The complete expressions for the matrix element as well as other intermediate steps 
in the calculations are given in Appendix A). This yields the following relativistic momentum space integral 
wave equation for stationary states of three identical particles: 

Fs{pi,2,3) (Wpi -I- + Wpj -E) = 

J M,2 3^3,3(p'l,2,3.Pl,2,3)^5(p'l,2,3) + j M,2,3 ^3,3 (pi,2,3> Pl,2,3) -^S (p'1,2,3) ■ (35) 
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Figure 2: The three-chion propagator vertex corresponding to the cubic interaction kernel £3,3 equation 
(37). The two propagators on the left should actually overlap (impossible to draw) such that they are 
perpendicular to the direction of time. 

The relativistic kernels of the interactions (i.e. the relativistic momentum-space inter-particle potentials) are 
given by 

3^3,3(p'l,2,3i Pl,2,3) 
<^3,3(p'l,2,3>Pl,2,3) 



where the summation is on the six permutations of the three indices 1,2 and 3. As is evident from the 
covariant factors on the RHS of equation (36), the kernel 3^3,3 corresponds to the three inter-particle one- 
chion exchange interactions. There arc no virtual annihilation terms since only particle operators (and no 
antiparticle operators) are present in the trial state (33). The kernel C3,3, equation (37), corresponds to the 
non-linear interaction term ^^/^ (details are given in Appendix A). It is similarly evident, from the covariant 
factors on the RHS of equation (37), that the Fcynman diagram corresponding to the cubic kernel C3^3 
contains a three-chion propagator vertex and is shown in Figure 2. Note that the three particle trial state 
does not probe the quartic interaction term 'H/g equation (16), i.e. {'9 3] Hj^ \^s) = 0. 

Equation (35) (with its kernels equations (36) and (37)) is one of our principal results. It is a relativistic 
wave equation for stationary states of a system of three identical scalar particles. It can describe purely 
bound states of the three particles or elastic scattering among them, but not processes that involve the 
emission or absorption of chions. The relativistic kinematics (i.e. kinetic energy) of the system is described 
without approximation, but the relativistic dynamics (i.e. potential energy) is described only at the level of 
one-chion exchange between the particle pairs by the relativistic kernel 3^3,3, along with a three-particle first 
iterative order correction C3.3 due to the non-linear interaction term T-Lj^ . 

The solution of the relativistic three particle equation (35) (even the non-relativistic limit, which we 
discuss below) is a challenging task which we shall not undertake in this work. No exact solutions are 
possible, therefore approximation methods (perturbative, variational, numerical or other) must be used. 
This shall be left for a future investigation. 

It is of interest and instructive to consider the non-relativistic limit of equation (35), in particular the 
coordinate representation of the interactions. Moreover, the non-relativistic description can be used for 
systems of heavy particles. 

The non-relativistic limit of equations (35)-(37), just as for the particle-antiparticle system, is obtained by 
assuming << and Fourier-transforming to coordinate space. From equation (36), the non-relativistic 
inter-particle potential for the three particle system due to the Yukawa interaction is, as expected, 

r g-A»a;i2 g-^a^is g-A»2:23 ^ 

Vy(xi,2,3,/^) = VY{Xij,f,) = -ag\ — + — + — \. (38) 

X12 Xi3 X23 



8(27r) 
8(27r)e 



7r)3 



^(Pi + P2 - Pn - Pt'J ^{P'a - P*3 



1 



»1,»2,23 
6 

E 



P.2 [m^ - (K -Pii)^ 

^(P'l + P2 + P-s - P»i - Pj2 - Pis ) 



(36) 



- {p'l +P2- Pn - Pt2Y - ip'i - PtiY - (P2 - Pi2Y 



(37) 
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Unlike in the particle-antiparticle case (32) , no virtual annihilation delta function contributions appear in this 
expression since we are dealing with a system of three identical particles only {i.e. there are no antiparticles). 
The non-relativistic cubic potential term that follows from equation (37) is 

,w , 3 /• , n-e""^'"'%i+q2 + q3) 

3 g^rj 

where a„ = — — — ^ — 5- is a coupling constant with dimensions of mass (see Appendix A for details). Vc can 
2567r°Tn'^ 

be simplified to a three-dimensional quadrature (see Appendix A for details): 

/g-/i|xi+x| g-/n|x2+x| g-/n|x3+x| 
|X1 + X| |X2 + X| |X3 + X| 

This is an overall well behaved convergent integral for any fi > 0. However, it cannot be evaluated analytically 
in general (at least we do not know how to do so). 

The expression for Vc, equation (40), with the substitution v = x + xi, can be written as 

/e-/*KI Q-f^W+^2i\ g-;i|v+x3i| 
|v| |v + X2l| |v + X3i| 

where X21 = X2 — xi and X31 = X3 — xi. As shown in Appendix C, we can also express Vc in the form 

Vc{^ij,H) = - ir'ar, / d/3i,2,3 ^ o o ^ o a ^ o o,s/2 ^^P ' ..o 0^0 o ^ o ' (42) 







(/3l/32 + /3l^3+/32^3)^/' "^V 4(A^2+/3i^3+/32^3)y 



which shows explicitly that Vc depends only on intcr-particlc distances Xij = Xji = |xj — Xj|. In light 
of equation (42), we conclude that the total potential energy function V = Vy + Vc depends only on the 
inter-particle distances, and is invariant under 3D rotations and translations of coordinates as we would 
expect of a closed system. Therefore, we shall also use the notation Vc{xij, ji) in what follows. In addition, 
equation (42) is suitable for numerical evaluations of Vc for arbitrary values of Xij . 

Unfortunately, it is not possible to plot inter-particle potentials of three independent variables X12, X23 
and xi3, even though they can be worked out numerically for arbitrary Xij. For this reason, we will plot one 
and two dimensional sections to get some idea of the shape of the cubic term Vc- 

We proceed to examine two particular cases for which Vc can be evaluated analytically. Then, we show 
surface plots of Vc for arbitrary xtj which we calculated numerically. 

Case 1: V(xi,2,2,m) = V{x2i,fi) 

The first case we examine is when /x > and X2 = X3 so that there is only one inter-particle distance, 

X21 = xi2 = |x2 — xij = |x3 — xi|, to be concerned with. The integral for Vc, equation (41), in this case is 
expressible in terms of the special function, the exponential integral, defined by 

roo t 

El (2) = J —dt, (43) 

namely 

Vc{^i,2,2,fi) = Vc{x2i,n) = -^^ie-^-''(ln(3)-Ei(a;2i/i)) +e^-''Ei(3a;2i/i) l (44) 

X2l/i [ ^ ^ J 

We choose to express length in units of the Bohr radius — ^— and, correspondingly, the energy in units of 

mag 

mai. Equation (44) is then written in terms of the dimensionless variables r = X2\Tna„ and M = ^ . 
y ^ ' " mag 

Suppressing the singularity at X2 = X3 in the Yukawa term (38), the total inter-particle potential V = Vy+Vc 

for the three particle trial state in the case when /U > and X2 = X3 is 

y(xi,2,2,/x>0) = y(r,M) 

f p-rM /p-rM p-r M rM \\ 

= »»^ -2— — 1„(3)-_E.(.M) + — E,(3rM) , (45) 
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47r**a 12r] 

where ki = 5^ = is a dimensionless constant. Figures (3) and (4) are plots of V{r,M)/ma^, 

ma^g g 

equation (45), as a function of r for the different values of M, with ki — and Ki = 0.1 respectively. It is 
apparent that the contribution of Vc lowers the overall potential with decreasing values of M. 

The behaviour of the inter-particle potential can be better understood by looking at the two extreme 
cases of the separation distance r. For small r, expanding equation (45), we obtain 

V(r,M)^-- + Ki\n(3rM) + CM, (46) 
r 

where Cm = 2 M + ki (7 — 1) is a finite constant with 7 = 0.5772 being the Eulcr's constant. The leading 
— 2/r term in the above expansion indicates that the inter-particle potential is dominated by the Yukawa 
term Vy for small r. To see the large r behaviour, we determine analytically the ratio Vc/Vy- 

^ = IM ( M3) - El (r M) + e^ El (3 rM)). (47) 

lnf3) Hi 

The large r limit of this ratio is the constant — —- — which implies that the cubic term Vn has the identical 

^ 4M 
large r dependence as the Yukawa term Vy. From this we see that the effect of Vc is to strengthen the 

attractive Yukawa-like potential, i.e. to increase the binding energies of bound states. 

Case 2: V{x2i = X31 = X23 = A, /x) = V{A, fi) 

Another case for which Vc can be expressed in closed from is when the coordinates are at the vertices 
of an equilateral triangle (see Appendix A for details). With the previous choice of the units of length and 
energy, we have 

V{A,ii) = V{r,M) 

= ""^"n^^^o (A/3. + A/33 + /3./33)3/^ ^"H"4(/3i/32 + /3i/33 + /3./33)j | ^''^ 

where r = Amag is the dimensionless inter-particle distance, A = |x2 — Xi| = |x3 — xi| = |x3 — X2I and 

M = — — is the dimensionless mass parameter of the mediating field. Here, we evaluated Vc numerically. 

mag 

The numerical integration was performed with the GNU Scientific Library. A plot of equation (48) is given 
in Figure 5. The curves show that V{r, M) approaches zero asymptotically from below for all M > 0. 

Ccise 3: Arbitreiry inter-particle distances 

Finally, we examine the inter-particle potential for arbitrary inter-particle distances. There are no ana- 
lytical solutions available for V^, equation (40), and the potential is calculated numerically with the help of 
the GNU Scientific Library. Using the parametrization given by equation (42), we obtain 



V{xij,iJL) = V{rij,M) = -ma 



2 



g-r23M Q—risM 



ri2 r-23 ri3 



^1 ' rz 

^4nJo (/3i/32+A/33 + /32/33)3/2 ^ \ W2 + l^iPs + W J j ^ ' 

where rij = Xij mag are the dimensionless inter-particle distances and Xij = |xi — Xj|. Four surface plots 
of equation (49) for diflterent values of ri3 are given in Figure 6. Note how the inter-particle potential rises 
with increasing separations among particles but does not cross the zero plane. 

We see that the cubic T-Lj^ interaction term of the Hamiltonian affects the inter-particle interaction 
potential V through the additional attractive term Vc. As a consequence, the energy spectrum for a three 
particle system will be lowered by the inclusion of the cubic term "H/^ . 
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10-^ 



10 



Figure 3: The inter-particle potential V{r, M) / ma^ of equation (45) as a function of r = 0:21 mag for 
M = 10-1° with Ki=0 (i.e. Vy only) and m = 0.1 (i.e. Vy + Vc). 




Figure 4: The inter-particle potential V{r, M) / ma^ of equation (45) as a function of r = 0:21 mag for 
Ki = 0.1 and various values of M. 



Four Particle State 

Since the throe particle trial state (33) does not probe the quartic Hj^ interaction term of the Hamiltonian 
{i.e. (^'ajifal^'a) = 0), wc consider a system of four identical particles and examine how both non-linear 
terms of the Hamiltonian, T-Lj^ and T-Li^ , affect the inter-particle potential. 

The four identical particle trail state analogous to equation (33) is given by 

1*4) = J dpi.A F(pi..4) at(pi) at(p2) at(p3) a\p4) |0). (50) 
For this four particle trial state to be an eigenvector of the momentum operator (25), we require that 
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Figure 5: The inter-particle potential V{r, M) / mag of equation (48) for ki = 0.1 and the indicated values 
of M = ^ , where r = Amag and A = X2i = X31 = X23- 



P 1^4) = Q |^'4) which can be achieved with the choice -F'(pi..4) = S{pi + P2 + Ps + P4 — Q) /(pi,2,3)- 
Thereupon, similar to the particle-antiparticle and the three particle cases, the wavefunctions will be of 
the form where the centre of mass motion is completely separable for this identical particle system. We 
shall write everything in terms of the completely symmetrized function Fs because the trial state (50) is 
completely symmetric under interchanges of the momentum variables. The symmetrized function is 

24 

Fs{pi.a) = 5Z PiPh,i2,i3,iJ^ (51) 

where the summation is on the 24 permutations of the indices 1, 2, 3 and 4. We calculate the matrix element 
(\I/4| H — E \ work out the variational derivative with respect to F* and set it to zero (see Appendix B 
for details). This leads to the following four identical particle relativistic equation for the function Fs (in 
momentum space): 

•Fs(pi..4) (Wpi +Wp2 +Wp3 +Wp4 -E) ^ j 2 y4,4(pi..4>Pl..4)-Fs(pi 4) 

+ j rfpi..4C4,4(pi..4,Pl..4)i^s(p'l..4)+ j dl?\„i e4,4(p'l..4,Pl..4)Fs(pi..4)- (52) 

The Yukawa and cubic interaction kernels 3^4,4 and C4,4 are similar in structure to those of the three particle 
trial state except there is dependence on an extra momentum coordinate (see Appendix B). The relativistic 
quartic interaction kernel is 

^ fi X 5^cr 1 5(pi +P2 + P3 + P4 - Pi - P2 - P.'i - P4) 

y4,4(.Pl..4)Pl..4j = -TT" /r, XQ 7, / 

^^^> il U V^Pi ^P2 '^PS '^Pl '^PS '^P4 



1\ ,12 ,13 ,14 

1 1 



- ip'i +P'2 +P3 -Pii -Pi2 -Pi^Y - (P'l -Pii) - {P'2-Pi2f - {P'3 -Pi^Y 



>, (53) 
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Figure 6: The inter-particle potential V{rij,M) / ma'^ of equation (49) for ki = 0.1, M = 
ri3 = 0.01(top left), 0.1(top right), 1.0(bottom left), 10(bottom right) where Vij = Aij mag. 
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Figure 7: The four-chion propagator vertex corresponding to the quartic interaction kernel Q4_4 equation 
(53). The two propagators on both sides should actually overlap (impossible to draw) such that they are 
perpendicular to the direction of time. 

where the summation is on the 24 permutations of the indices 1,2,3 and 4. This contribution to the potential 
corresponds to a four-chion propagator vertex shown in Figure 7. 

It is not our intent to solve the four-particle rclativistic equation (52) in this paper. No exact solutions 
are possible, just as with the three-particle rclativistic equation (35). Instead, we proceed to examine the 
non-relativistic limit of the interactions in the coordinate representation. 

In the non-relativistic limit, the Yukawa and cubic kernel reduce just as their three particle trial state 
counter-parts. There are six Yukawa interaction terms which can be worked out in analytical form: 



VY{-X-i.A,fJ-) = -ag 



Xl2 



Xl3 



Xl4 



X23 



-fJ.X2i 



X24 



-11X34, 



Xzi 



(54) 



where aa 



is the dimensionless coupling constant and Xi 



the inter-particle 



distances as before. 

The cubic interaction kernel for the non-relativistic four particle case reduces to a sum of four terms, 
i.e. one for every three-way interaction: 



Vc(xi..4,m) = -"r, / (iqi.. 



+ 



*\ (M^ + qD(M^ + qi)(M^ + qi) 

n^c-'q.:--,5(qi+q3 + q4),5(q2) 
(/x2 + q2)(^2+q2) (^2 + q2^ 



, 'q'-'J(qi + q2 + q4)'^(q3) 
(/z2+q2)(^2+q2)(^2 + q2) 

n-e-'q.--J(q2 + q3 + q4)%l) ] .... 
(M2+qi)(/.2 + q|)(/,2 + q2) J- l^^i 



where a„ 



^rj — K — T is a coupling constant with dimensions of mass (see Appendix B for details). The delta 

function with a single momentum variable in every term indicates that the particle carrying that momentum 
is a "spectator" of the three-way interaction. To obtain this expression we followed similar steps as those 
leading to equation (39). 

The quartic interaction kernel for the four particle system in the non-relativistic limit reduces to 



'5(qi + q2 + qa + q4) 



(M^ + qf) (M^ + qi)(M^ + ql)(M^ + ql) 



2\' 



(56) 



where a„ 



3gV 



is a dimensionless coupling constant. To obtain this expression wc used similar steps 



1024 TT^m^ 

as those leading to equation (55). The cubic Vc and quartic Vq contributions to the inter- particle potential 
are both symmetric under 3D rotations and translations even though it is not readily evident from equations 
(55) and (56). In Appendix C, equations (96) and (99), we show explicitly that they indeed depend only on 
the inter-particle distances, i.e. Vc = Vc{xij) and Vq = VQ{xij). 
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The expressions for Vc and Vq can be simplified to three-dimensional quadratures: 

g-M|xi+x| g-Allx2+x| g-/i|x3+x| g-/j|xi+x| g-;ilx2+x| fJ.\x4,+x\ 



1x1+ x| 1x2 +x| IX3-I-XI 1x1+ x| 1x2+ x| 1x4+ X| 

g-/i|xi+x| g-/^|x3+x| g-/i|x4+x| g-/^|x2+x| g-/^|x3+x| g-/i|x4+X 



|X1+X| |X3+X| |X4+X| |X2+X X3+X X4 ■ 



(57) 



/g-/x|xi+x| g-/i|x2+x| g-;i|x3+x| g-/j|x4+x| 
dx -; 1-- 1--; 1--; 1- . (58) 
|Xi+x| |X2+X| |X3+X| |X4+X| ^ ' 

Here, analogously to the three particle case (40), the integral expression for Vc, equation (57), is convergent 
for /i > 0. The integrand in Vq, equation (58), goes as |x|~^ for large |x| and the integral remains finite for 
/i > 0. As wc shall see below, the contribution of Vq to the total intcr-particlc potential of the four identical 
particle system is particularly significant for small separations, depending on the values of the coupling 
constants. 

The expressions for Vc and Vq, equations (57) and (58), analogously to equation (41), can be written as 

g-A»kl I g— /'|v+X2l| g— /:i|v+X3i| g— /n|v+X2l| g-ju|v+X4i| 

+ 



/e-A«|v| ( Q-i 



+ X2l| |v + X3i| |v + X2l| |v-|-X4i| 

g-/^|v+X3i| g-/j|v+X4i| g-A'|v+X32| (,-;<|v + Xl2 



V + Xsi V-I-X41 V-I-X32 V-I-X42 



(59) 



/e-A»KI e-A'k+x2i| g-Ai|v+x3i| g-/j|v+x4i| 
'^'^ ~n — \ — \ n — ] n — \ T' (^o) 
|v| |v + X2i| |v-|-X3i| |v + X4i| 

where Xjj = Xj — Xj . 

It is not possible to obtain analytical expressions for Vc and Vq in general and one must resort to 
numerical evaluation of these integrals. However, as for the three particle system, there is a particular solvable 
case, namely when xi = X3 and X2 = X4, which shows the general features of the total potential energy. 
With this restriction on the coordinates there is only one inter-particle distance 3:21 = |x2 — xi| = |x4 — X3I 
in the problem. 

Special Case: V(xi,2,i,2, M > 0) = V{x2i,n > 0). 

For the case when /x > the cubic potential term results in four copies of what was previously obtained 
for the three particle case, namely, 

Vc(xi,2,i,2,M > 0) = Vc{x2i,n> 0) = Je— " [ln(3) -Ei (a;2iM)] + e"^^'* Ei (3a;2iM) }■ (61) 

X21IJ' I J 

The quartic potential term Vq(xi^2,i,2) M > 0), equation (60), can be reduced to an evaluation of a single 
quadrature (see Appendix B for details): 

Vq(xi,2,i,2,M > 0) = VQ{x2uf^> 0) =?^J r" — e-2^'^|Ei(2M(a;2i -t;)) -Ei(2/x(a;2i +t;))| 

X21 [ Jo ^ L J 

+ ^^~^''"{MMv-X2i))-MMx21+v))^^. (62) 

The integrals in the expression above can not be expressed in terms of common analytic functions, thus we 
evaluated them numerically using the GNU Scientific Library. 

The total potential energy V = Vy + Vc + Vq in the units of ma^ as a function of r = a;2i mag for the 
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case when /U > 0, suppressing two Yukawa singularities, is 
l^(xi,2,i,2,/U > 0) =F(a;i2,M>0) 




El (2(rM - w)) - El (2(rM + w)) | 

^ V ^^^"^ ~ ""^^^ ~ ^^^""^ ^ ""^^ }) } ' ^^^^ 

where w = t; is a dimensionless integration variable and M = ^ is the dimensionless mass parameter 

mag 

as before. Two plots of equation (63) arc shown in Figure 8 for the parameter values ki = 0.01 and 
K2 = 0.1, 0.005. The choice K2 = 0.1 yields a potential which is repulsive for small separations and possesses 
a trough. The trough's depth increases with increasing values of the mediating field mass parameter M. For 
large separations and both values of K2, the potential approaches zero from below without crossing the zero 
of energy. 



Concluding Remarks 

The goal of this work was to derive and analyze relativistic wave equations, including the interactions 
kernels, for systems of three and four particles in theories that contain non-linear mediating fields. We 
investigated a non-linear theory (1) in which scalar "matter fields" (f> and cj)* interact via a real non- linear 
scalar field x of the Higgs type. We reformulated the theory by using the formal solution (4) of the classical 
equations of motion for the field x- This enabled us to express the interaction terms of the Hamiltonian in 
terms of the x-ficld propagator and the matter fields. Since the integral equation (4) cannot be explicitly 
solved for x, we used an iterative approximation and obtained the Hamiltonian density (12), which does not 
contain x explicitly but only implicitly through the X-^eld propagator. 

We used the variational method to derive approximate few-particle eigenvalue equations for the quantized 
version of the theory. To this end we used Fock-space trial states, equations (33) and (50), to calculate the 
matrix elements (^"3! H l'^^) and (^4] H 1*4) for three and four identical-particle systems. Thence, we used 
the variational principle (19) to derive the relativistic wave equations (35) and (52) for stationary states of 
the three and four identical particle systems. The kernels of these integral equations represent the relativistic 
inter-particle interactions in momentum space. 

In order to get a picture of the inter-particle interactions in terms of familiar concepts, we considered the 
non-relativistic limit of the few-particle equations (Fourier transformed to co-ordinate space). 

The particle-antiparticle system described by the trial state (27) served as a preamble. For this system, 
the contributions of the non-linear, cubic and quartic, mcdiating-field terms rj x^ and a x'^ in equation (1) 
are not sampled by the trial state (27), and so the non-relativistic particle-antiparticle potential is just 
the Yukawa plus virtual-annihilation contact potential (32). To probe the eff'ect of the non- linear terms, 
one would need to consider a more elaborate trial state than (27), such as one containing more Fock-space 
components, for example 

|*2> = |*2> + 1*4) 

= J rfpi,2 i=^(pi,2) at(pi) 6t(p2)|0) + J dpi.A G(pi..4) at(pi) bHp2) a^ps) 6+(p4)|0). (64) 

The trial state (64) leads to coupled multidimensional integral equations for the coefficient functions F 
and G. These equations would need to be solved approximately (say, variationally), which is a tedious, if 
straightforward task. The extraction of information about the particle-antiparticle potential would be an 
additional challenge. 

Reformulation of the theory enabled us to use simple three and four particle trial states, equations (33) 

and (50), to derive relativistic wave equations and interaction kernels which include novel additional terms 
due to the non- linear terms of the Hamiltonian density and T-Lj^ . The non-relativistic limits of these 
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Figure 8: The inter-particle potential V{x2i,H > 0) of equation (63) in units of ma^ as a function of 
r = X21 mag for the indicated values of M, m = 0.01, and K2 = 0.1 (top) and 0.005 (bottom). 



interactions, in coordinate space, are three and four-point inter-particle potentials which are shown to depend 
on the inter-particle distances only. 

Unfortunately, the potentials (40), (59) and (60) due to the non-linear interaction terms cannot be eval- 
uated analytically for arbitrary inter-particle distances, and hence must be evaluated numerically. However, 
analytic expressions are obtainable for particular, restricted cases, for arbitrary values of the mass /x of the 
mediating field quantum, and these provide a good picture of the general form of the potentials. 

For the three particle system, the analytic expressions of the inter-particle potentials with certain coor- 
dinate restrictions are given in equations (45) and (48) with > 0. We find that the Yukawa potential is 
augmented by the cubic term Vc as illustrated in Figures 3, 4 and 5. The potential does not rise above the 
reference line of the rest mass {i.e. 2m) as the inter- particle separation increases but simply approaches zero. 
This indicates that the potential does not support bound states with energies greater than the rest mass of 
the constituents. 

For the three particle system, the numerical evaluation of the potential in the case of arbitrary inter- 
particle distances is performed where one of the three inter-particle distances is kept constant. These surface 
plots are displayed in Figure 6. Similarly, they show that the inter-particle potential does not support bound 
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states with energies larger than the sum of the masses of the constituent particles. 

For the four particle trial state with n > 0, the analytic expression of the potential for the particular case 
Xi = X3 and X2 = X4 is given in equation (63). The quartic term can alter the behaviour of the potential 
for small values of the inter-particle separation depending on the coupling constants, as is shown in Figure 
8. High values of the mass parameter M are seen to be responsible for the development of a trough in the 
potential, provided that the coupling constants ki and K2 are such that the potential is repulsive for small 
r. Unfortunately, it is not possible to give graphical representations of the potentials for arbitrary values of 
the coordinates. 

In this work we have concentrated on the nature of the intcr-particlc interactions in a scalar model with 
a non-linear Higgs-type mediating field. The non-relativistic limit was considered in some detail. We shall 
leave for the future the challenging task of solving the three and four particle relativistic and non-relativistic 
equations derived in this paper, as well as the analysis of the particle-antiparticle system with more elaborate 
trial states. 



Appendix A: Three Particle State 

In this section we present some intermediate steps of the derivations with the three particle trial state 
(33). All matrix elements were calculated using Maple. 

Derivation of the kernels 

The matrix element in the Schrodinger picture for the three identical particle trial state (33) is given by 



(*3| H-E 1*3) = (*3| + Hi, +Hi,-E 1*3) 



where the contributions are 



{^3\H^-E\^3) = j dpi,2,3(wpi + UJp, + U)p, - E) F*{p[^2,3)Ps{Pl,2,3), 

'^Pl,2.3^Pl,2,3 



(*3|^?/J*3) 



8 (27r)3 J ^ujp'^ujp'^ujp^ujp. 
X 5{p[ + P2 - pi - P2) ^(Ps - P3 



^S(P'1,2,3)^^S(P1,2,3) 
1 



(M-ali?/, 1*3)=- 



8(27r 



- ijp'i -piY 



(65) 



(66) 



(67) 



<^P'l,2,3«^Pl,2,3 



X (5(pi + P2 + P3 - Pi - P2 - P3) 

1 1 



^^S(P'1,2,3)^^S(P1,2,3) 



1 



.M^ - {p'l +P2 - -Pi -P2Y - (p'l -PiY IJ-"^ - {p'2 -P2Y 
In working out the variational derivative the following identity was used: 

6 



1,2,3; 



^-F(qi,2,3) 



X] ^(Pii - qi) <5(Pi2 - q2) <5(pi3 - Qs), 



(68) 



(69) 



where the summation is on the six permutation of the indices ii. 12 and i^. 

The relativistic kernels for the three particle trial state are given in equations (36) and (37) in the body 
of the paper. In the non-relativistic limit the kernels become 



^3,3(p'l,2,3>Pl,2,3) = 



8 (27r)3m2 

(5(pi-|-p^-pi-p3)(5(p^-p2) , ^(P2 + P3 - P2 - P3) ^(P'l - Pi) 



'^(Pl + P2 - Pi - P2) '?(P3 - P3) 
+ (Pi - Pl)2 



+ 



^3,3(p'l,2,3>Pl,2,3) 



M' + (P^-P3)2 M2 + (P^-P2)2 

3 g'^v ^(p; + p[, + Tp'^ - Pi - P2 - P3) 



4 {2Tl fw? (/X2 + (pi - pi)2)(/x2 + (p^ - P2)2)(m2 + (Pjs - P3)2) 



(70) 
(71) 
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where we used the symmetry property of Fg to simphfy the expression. 

3 

To obtain equation (39) we multiplied equation (71) by JJ^ e'^''^' and shifted the variable of integration 

i 

p[ —Pi = Qi- Integrating over the vectors leads to equation (40). This integration is performed using the 
standard technique where the radial integral is evaluated in the complex plane using Cauchy's integration 
formula. 

Derivation of the cubic potential energy Vc for X2 = X3 

The cubic inter-particle potential for the three particle trial state (33) in the case when /i > and X2 = X3 
is determined from the following expression: 

/g-M" g-2M|v+X2i| 
, — : 12- ' ('''2) 
V |v + X2l|^ 

where X21 = X2 — Xi and the non-vector notation is understood to mean the magnitude of the vector. After 
the trivial azimuthal integration we obtain: 

Vc(a;2i,M) = -27r% / dvv^ dw — . — -, (73) 

Jo J-i ^ |v + X2l| 

where w = cos 9 with 9 being the polar angle of the vector v. The polar integration is performed using the 
substitution g'^ = + i;^ + 2.r2i w w, with qx — {x2\ + v) and Q2 = |.'£'2i — v\ as the upper and the lower 
limits of integration. Thereupon, the integral of equation (73) can be written as 

Vc{x2uf^) = r^e-^i^^ r dge-^^. (74) 

Jo V jg^ 

We integrate over the variable g, split the interval of integration accordingly and end up with the following 
radial integral: 

Vc{x2i,n) = Jc-'^-- f — (e-'"' - e-3'^^') + (e"-- - e"''^-) H — e-^'^A. (75) 

1^x21 [ Jo V ' J^^^ V J 

The integral is expressible in terms of the exponential integral defined by equation (43). Thus, the cubic 
potential term for the three particle trial state in the case when fj. > and X2 = X3 evaluates to 

^"c(a;2i,/^) = -^^|e-"-qin(3)-Ei(x2i/x)]+e^-''Ei(3:r2iM) \, (76) 

X21l^ [ J 

which is equation (44). 

Appendix B: Four Particle State 

In this section we present some intermediate steps of the derivations with the four particle trial state 
(50). All matrix elements were calculated using Maple. 

Derivation of the kernels 

The matrix element in the Schrodinger picture for the four identical particle trial state (50) is given by 
(*4| H-E 1*4) = (*4| + Hi, + Hi, + Hi, - E 1*4), (77) 
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where the contributions are 



(*4| H4,-E 1*4) = y dpi.A F*{pi.a) Fs(pi..4) (wp^ + cop^ + + - E) , 

0'p'l..4''Pl..4 



X ^(pi + P2 - Pl - P2) ^(Ps - Ps) 5(P4 - P4) 



and 



8(27r)6 J ^Wp'^Wp^Wp^Wp^Wp.Wpg 

X ^(pi + p'2 + Ps - Pl - P2 - Pa) ^(P4 - P4) 

1 1 1 

- {p'l +P'i-Pi- P2Y - {P'l - PiY - (P2 - Pa)^ 

dp'i 4dpi..4 



(vl/4|if/3|M/4) = 



16(27r)9 y ^Wp.Wp/Wp/Wp/w PjWpjWpgWp^ 
X 5(p'i + P2 + Ps + P4 - Pl - P2 - P3 - P4) 



F^(p'l..4)F5(pi..4) 



- (p'l + P'2 + Ps - Pl - P2 - Ps)^ 
1 1 



1 



.M^ - (p'l - Pl) - (P2 - P2)^ - (Ps - Ps)^ 
The following identity was used in working out the variational derivative: 



'^i^s(pi..4) 

^F(qi..4) 



= ^(Pii - qi) '5(Pi2 - qa) ^(Pi3 - qa) ^iPu - 



where the summation is on the 24 permutation of the indices ii, 12, is and 14. 

The relativistic Yukawa and cubic interaction kernels for the the four particle trial state (50) are 



3^4,4(p'i..4,Pl..4) = - 



16(27r)3 

24 



E 



^(P'l + P2 - P^i - Pm) 'Kp's - Pia) '^(P4 - V>ii) 



WpiWp/Wp,^Wp,^ 



- (p'l -PiiY 



9^V ^(Pi +P2 + P'i - Pa - P'2 - P' J '^(P4 - Pu) 



C4,4(p'i..4)Pl..4) 



8(27r)6 



E 



1 



WpiWp^Wp^Wp,^Wp,^Wp,3 

1 1 



.M^ - (p'l +P'2 -P«i -P«2)^ - (p'l -Pn)^ - (p'2 -Pi2)^ 
The relativistic quartic kernel for the four particle trial state is given in equation (53). 
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In the non-relativistic limit the interaction kernels for the four particle trial state reduce to 

^2 



5^(pi..4>Pi..4) 



4(27r)3m2 







+ (pi-Pl)2 




+ 




+ (Pi - Pi) 


2 




-^(Pi+P^ 


- Pi - 


-P4)'5(P2 -P2)^(P3 


-Ps) 


+ 


^(P2 + Ps 


-P2 -P3)'5(pi - 


Pi)'5(p1 


-P4) 






+ (P'l - Pl)2 






+ (P2 - P2 






<5(P2+P4 


- P2 


- P4)<5(pi - Pl)<5(P3 


- Pa) 


+ 


<5(P3 + P4 


- P3 - P4)<5(pi - 


-Pi)(5(p^ 


- P2) 



M' + (P^-P3)2 



C(p'l..4>Pl..4) 



3 g'^ri 



4(27r)6TO3 



+ 



+ 



'^(P'l + P2 + P3 " Pi - P2 - P3) '^(P4 - P4) 
(^2 + (p'l - Pl)2)(/Z2 + (p^ - P2)2)(^2 + (p, _ p3)2) 

^(P'l + P2 + P4 - Pi - P2 - P4) ^(Ps - Ps) 



(/Z2 


+ (pi-Pi)2)(M2- 


Kp'2 


-P2)2)(/i2 + (p^-p4)2) 




'^(P'l+P3 + P4- 


Pi 


P3 - P4)<5(P2 - P2) 




+ (p'i-Pi)2)(m2- 


Kp!i 


-P3)2)(/X2 + (P;-P4)2) 




^(P2 + P3 + P4 - 


P2 - 


P3 -P4)(5(p'i -Pl) 



(M^ + (P'2 - P2)2)(/X2 + (p^ - P3)2)(^2 + (p. _ p^)2) 



and 



<3(p'l..4)Pl..4) 



3gV 

'2(27r)W 



<^(P'1 + P2 + P3 + P4 - Pl - P2 - P3 - P4) 



(84) 



(85) 



(86) 



.(M^ + (p'l - Pl)^)(M^ + (P^ - P2)2)(m2 + (p(, - P3)2)(^2 + (p/ _ p^)2) 

Derivation of the queirtic potential energy Vq for Xi = X3 and X2 = X4 

The quartic inter-particlc potential term for the four particle trial state (50) in the case when /x > 0, 
xi = X3 and X2 = X4 is determined from the following expression: 



^q(xi,2,i,2,M > 0) = VQ(a;2i,A« > 



0) = a^TT^ j 



dv 



-2At|v| g-2Ai|v+X2i| 
|v + X2l|2 ■ 



Following similar steps as those in equations (72)-(74) leads to 



Vq(xi,2,1,2,M > 0) = VQ{x21,fi > 0) 



27r^cr f dv 



X21 



]^ -2tJ.v / l!t:p-2/ie 



dg 



(87) 



(88) 



where gi, g2 and X21 are given below equation (73). The result of the integration over the variable g can be 
expressed in terms of the exponential integral (43) and is given in equation (62). The remaining integrals 
have to be evaluated numerically. 



Appendix C: Gaussian Parametrization 

In this section we present an alternative method for evaluating the cubic and the quartic potential 
contributions Vc and Vq for the three and four identical particle cases. The invariance of the potentials Vc 
and Vq under rotations and translations of the coordinates becomes explicit in this method. 

We shall write the denominators of equations (39), (55) and (56) using a technique which enables us to 
perform Gaussian integration over the momentum variables. We use the following identity and the Gaussian 
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integration formula: 



d/3e-^'^=- for A>0, 
A 



/■OO 
dxi exp 
-OO 







' K' 


2 


LjXj - W 


= Det 


2^ 


exp 



-Kr/L,Lj-W 



(89) 
(90) 



where K is an invertible and symmetric matrix, L is a vector and is a constant. 

Three Particle State 

The three particle trial state (33) yields the cubic interaction kernel (39). After integrating over the delta 
function in (39) there remain only two momentum integrations. Using the identity of equation (89) on each 
factor separately, we obtain the following expression: 



Vc(xij,/i > 0) = -TT^ttr, y rf/3i,2,3 J dqi,2 e" 



(/i^ ) 01 - {n^+ql ) /32 - (/i^ + (qi +q2 ) fop-iqi -^21 p-iq2 -xsi 



(91) 



where the inter-particle vectors. Upon expanding the squares in the exponentials and 

defining q = qi + qi we identify a six dimensional Gaussian integral. The matrix K and its inverse K~^, 
the vectors and W for the Gaussian integral are given in block diagonal form 



K 



2(/3i+/33) 2/33 
2/33 2(/32 + /33 

L =i [ X21, X31 ] , 



1 



4(/3i/32 + ■ 
= (/3i+/32 + /33)/z2. 



-P2P3) 



2(/32 + /33) -2/33 
-2/33 2(/3i+/33) 



Applying equation (90) we arrive at 



Vc(Xij,A* > 0) = -TT'^a^ / d/3i,2,3 



-,i=(/3i+;82+/33) 



(A^2 + /3l^3 + /32/33)3/2 



exp 



/3ixii + ^2xii + /33xi2 

4 (/3i/32 + ^1/33 + /32/33) 



(92) 
(93) 

(94) 



This expression shows explicitly that the cubic potential term Vc depends only on the inter-particle distances. 
Numerical integration over the parameters /3i is required to complete the calculation. 



Four Particle State 

The four particle trial state (50) yields the cubic interaction kernel (55). The calculation of Vc for the 

four particle trial state is identical to that of the three particle trial state. There are basically four copies of 
the three particle results with different inter-particle distances involved. The matrix K and its inverse 
are identical for all four terms and the same as for the three particle trial state equation (92). The vectors 
Li pertaining to each term in equation (55) are 



Li=i[ X21, X31 ] , L2 = i [ X21, X41 ] , is = i [ X31, X41 ] , L4 = i [ X32, X42 



(95) 



where the subscript indicates the corresponding term in equation (55). The constant W is the same as in 
equation (93). Application of the Gaussian integration formula (90) leads to the result 



Vc(xy ^ > 0) = - a,,TT^ 



exp 



exp 



(/3i^2 + /3i/33 + /32/33)3/2 

/3ixf2 + fex^3 + /33XI3 

4(/3i/J2 + 13103 + 13203) 

^ixfa + 02^14 + ^3X^4 
" 4(^1^2 + /3l/33 + /32^3) 



exp 



exp 



/3lX?2 



- /32X?4 



03^14 



4(/3i/32 + 0103 

01X^3 + (^2^24 



■ 0203) 
- /33X34 



4(/3i/32 + (3i(33 



Prom equation (96) it is evident that Vc depends on the inter-particle distances Xij 



■0203), 

|xj - Xj I only. 



(96) 
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The four particle trial state yields the quartic interaction kernel equation (56). The calculation follows 
the same steps as for the cubic interaction kernel. Applying the identity of equation (89) leads to a 9 
dimensional Gaussian integral. The matrix K and its inverse 'm block diagonal form are 



K = 



2(/3i + ^4) 2/34 2^4 

2/34 2(/?2 + Bi) 2/34 
2/34 2/34 2(/33 + Bi) 



2p 



1234 



^234 -Mi -Ma 

-Ma /3i34 -Ma 
-Ma -~Ma Bi2a 



(97) 



where ^1234 = M2fi3-\- IBiMa-^ M2Ba-\- MzIBa and jBijk = (Bij + Bik-^ IBjk- The vector L and the constant 
W of the Gaussian integration are 

L = i(x4i, X42, X43), W ={pi+ p2+ Pi + Pa)p?. (98) 
Applying the Gaussian integration formula (90) and after some algebra we end up with the expression 



poo 

Vqixij,!! > 0) = a^TT^/^ / rf/3i..4 

Jo 

X exp 



(^234)3/2 

Ma^12 + MA^la + M3^1a + Ma^I^ + Ms^Ia + M2XIA 



,2 
■■14 

4/3i234 



(99) 



This expression shows explicitly that the quartic potential term Vq depends only on the inter-particle 
distances. Numerical integrations over the parameters ft are required to complete the calculation of Vc and 
Vq. 
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